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programming  models  with  optimization  problems  in  the  constraints 
have  been  discussed  in  a  previous  paper  [1].  A  computer  program 
for  solving  mathematical  programming  models  with  nonlinear 
programs  in  the  constraints  has  been  presented  in  a  subsequent 
paper  [2].  A  procedure  for  transforming  mathematical  programs 
with  two-sided  optimization  problems  in  the  constraints  into 
mathematical  programs  with  nonlinear  programs  in  the  constraints, 
enabling  solution  by  the  computer  program  of  [2],  has  been 
given  in  [3] . 

The  present  paper  formulates  models  of  defense  problems 
which  are  convex  programs  having  the  mathematical  properties 
treated  in  the  previous  papers.  The  models  include  several 
strategic  forces  planning  models  and  two  general  purpose  forces 
planning  models. 


DD  ,'.r.,1473  UNCLASSIFIED 

Security^  Cl  a  m  si  fi  cation 


UNCLASSIFIED 

Security  Clarification 


1  4 

key  words 

LINK  A 

LINK  B 

LINK  C 

ROLE 

WT 

ROLE 

*T 

ROLE 

W  T 

Optimization 

Defense  Models 

Strategic  Forces  Planning 

General-Purpose  Forces  Planning 

Two-Sided  Optimization 

Nonlinear  Programming 

Game  Theory 

Max-Min 

< 

I 

UNCLASSIFIED 

Security  Claaaification 


ABSTRACT 


Theory,  computation,  and  an  example  of  mathematical  programming 
models  with  optimization  problems  in  the  constraints  have  been  dis¬ 
cussed  in  a  previous  paper  [1].  A  computer  program  for  solving 
mathematical  programming  models  with  nonlinear  programs  in  the 
constraints  has  been  presented  in  a  subsequent  paper  [2].  A 
procedure  for  transforming  mathematical  programs  with  two-sided 
optimization  problems  in  the  constraints  into  mathematical  programs 
with  nonlinear  programs  in  the  constraints,  enabling  solution  by 
the  computer  program  of  [2],  has  been  given  in  f  3] . 

The  present  paper  formulates  models  of  defense  problems  which 
are  convex  programs  having  the  mathematical  properties  treated  in 
the  previous  papers.  The  models  include  several  strategic  forces 
planning  models  and  two  general  purpose  forces  planning  models. 


I .  INTRODUCTION 


Reference  [1]  presents  theory,  interpretations  and  an  example  of 

mathematical  programs  with  optimization  problems  in  the  constraints . 

The  first  result  there  deals  with  Problem (A):  find  vectors 

x  =  (Xp  xn)  and  v1  =  (v^,  v^  ),  for  i  =  1,  m,  to1 

i 

minimize  f(x) 
x  e  X 

subject  to 

h.(x)  =  min  g1(x,v1)  >0  ,  i  =  1,  m 

i 

v  eV 

If  f(x)  is  a  convex  function  of  x  on  a  convex  set  X  and  g1(x,  v1)  is 
concave  in  ’ x  on  X  for  every  v1  e  v1  and  for  i  =  1,  m,  then  the 

mathematical  program  is  a  convex  program.  In  order  to  computationally 
solve  this  problem  it  is  desirable  that  g1(x,  v1)  be  convex  in  v1  and 
the  set  V1  be  convex.  In  Reference  [1]  an  example  is  presented  and 

solved  in  which  g1(x,  v1)  is  concave  in  x  ,  and  min  g1(x,  v1)  is  a 

i  r,i 
v  eV 

linear  programming  problem  in  v1 . 

The  second  result  in  Reference  [1]  deals  with  Problem  (B):  find 
vectors  x  =  (x1?  ...,  xn),  v1  =  (v^,  ...,  v^),  and  u1  =  (uj,  ...,  u^  ), 

for  i  =  1,  . . . ,  m,  to 


1.  All  indicated  maximums  and  minimums  are  assumed  to  exist. 
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— 


— 


minimize  f(x) 
x  e  X 

subject  to 

^(x)  =  max  min.  g1(xJ  u1,  v1)  >  0  ,  i  =  1,  . . . ,  m  . 

u1ell1(x)  v1eV1 

If  f(x)  is  a  convex  function  of  x  on  a  convex  set  X  ,  g1(x,  u1,  v1) 
is  concave  in  (x,  u1)  for  every  v1  e  V1, 'arid  U1(x)  is  a  convex  set  for 
x  e  X,  for  i  =  1,  . ..,  m,  then  the  mathematical  program  is  a  convex 
program.  In  order  to  computationally  solve  the  problem  it  is  also 
desirable  that  g1(x,  u1,  v1)  be  convex  in  v1  and  the  set  V1  be 
convex. 

Reference  [2]  presents  a  computer  program  for  solving  Problem  (A) 

where  the  constraints  min  g  (x,  v  )  >  0,  i  =  1,  . ..,  m  contain  non- 

i  „i 
v  eV 

linear  programs . 

Reference  [3]  shows  how  Problem  (B)  can  be  solved  by  the  computer 
program  of  Reference  f*  2 ]  . 

The  purpose  of  this  paper  is  to  present  optimization  models  having 
the  structure  of  the  above  mathematical  programs  for  some  defense 
planning  problems .  Each  of  the  problems  that  follows  is  formulated  in 
terms  of  its  major  planning  variables.  Specific  functional  forms  are 
assumed  which  have  the  general  properties  possessed  by  the  physical 
process  being  modeled.  These  functional  forms  have  parameters,  the 
values  of  which  could  be  determined  from  empirical  data  or  detailed 
quantitative  descriptions  of  the  physical  processes.  This  paper 
concentrates  on  the  problem  formulation.  All  of  the  models  considered 
are  convex  programs. 
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II.  GENERAL  CHARACTERISTICS  OF  THE  MODELS 

Consider  Problem  (A).  The  constraint  set  for  the  mathematical 
program  is  equivalent  to 

{  x  e  X:  g1(x,  v1)  >  0  for  all  v1  e  V1  ,  i  =  1,  . . . ,  m  }  . 

This  formulation  renders  a  ready  interpretation  of  the  problem.  The  "outside 
optimizer"  chooses  forces  x  such  that  each  of  the ...  .m  effectiveness 
requirements  in  the  constraint  set  can  be  met  at  a  minimal  value  of 
the  objective  function,  no  matter  what  feasible  choice  of  v  ,  ...,  vm 
the  "inside  optimizer"  makes.  An  alternative  interpretation  is  that 
of  the  outside  optimizer  choosing  x  ,  followed  by  (or  simultaneous 
with)  the  inside  optimizer’s  choice  to  minimize  the  effect  of  x  . 

The  first  military  model  having  the  structure  of  Problem  (A)  is 
that  of  an  outside  optimizer  choosing  forces  x  with  minimum  cost  f(x) 
which  can  achieve  effectiveness  objectives  r.. ,  ...,  t  in  the  face 
of  optimal  allocations  of  v1  e  V  (i  =  1,  . . . ,  m)  by  one  inside  optimizer 
in  an  attempt  to  minimize  the  effects.  The  model  is  to  choose  x  and 


minimize  f(x) 

subject  to 

min  g1(x,  v1)  >  r^  ,  i  =  1,  . . . ,  m 
v1eV 

The  outside  optimizer  can  achieve  his  objectives  against  any  of  the 
various  uses  of  one  inside  optimizer’s  resources  V  . 
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The  second  military  model  having  the  structure  of  Problem  (A) 
is  that  of  an  outside  optimizer  choosing  forces  x  with  minimum  cost  f(x) 
which  can  achieve  any  one  of  the  objectives  r^,  . ..,  in  the  face 
of  optimal  allocations  of  v1  e  V1  (i  =  1,  m)  by  m  inside 

optimizers.  The  model  is  to  choose  x  and  v  ,  . ..,  vm  to 

minimize  f(x) 


subject  to 

min  g1(x,  v1)^-  r.  ,  i  =  1,  m 

i  „i 
v  eV 

The  outside  optimizer  chooses  x,  and  m  inside  optimizers  choose 
v1,  .  ..,  vm.  The  outside  optimizer  will  be  capable  of  achieving  his 
objective  in  any  one  conflict  in  which  he  engages.  One  example  might 
be  the  outside  optimizer  planning  to  fight  and  win  any  one  conventional 
war  but  not  more  than  one  simultaneously. 

The  third  military  model  having  the  structure  of  Problem  (A)  is 
that  of  an  outside  optimizer  choosing  forces  x\  ...,  xm  with  minimum 
cost  f(x\  ....  xm)  which  can  achieve  all  of  the  objectives  rn ,  ....  r 

in  the  face  of  optimal  allocations  of  v1eV1  (i  =  1,  ...,  m).  The  model 

,  1  m  ,  1  m 

is  to  choose  x  ,  . . . ,  x  and  v  ,  . . . ,  v  to 

minimize  f(x  ,  ...,  x  ) 


subject  to 


min 

v1eV 


i 


i,  i 

g  (x  . 


V 


)  >  ri 


i  =  1, 


m 
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There  are  two  interpretations  of  this  problem.  There  may  be  one  inside 
optimizer  who  divides  his  resources  into  V  ,  V™,  or  there  may  be 

m  inside  optimizers  with  resources  v  ,  . ..,  V™.  In  either  case,  the 
outside  optimizer  can  simultaneously  achieve  all  of  the  objectives. 

Consider  Problem  (B).  The  constraint  set  for  the  mathematical 
program  is  equivalent  to 


|  x  e  X:  min  g1(x,u1,v1)  >  0  for  some  u1  e  U1(x)  ,  i  =  1,  m}  . 

i  „i 
v  eV 

Again,  this  formulation  provides  a  ready  military  interpretation.  The 
outside  optimizer  chooses  forces  x  .  For  this  choice,  there  is 
a  feasible  set  of  alternative  uses  of  these  forces,  in  achieving  an 
effectiveness  objective  i  ,  namely  U1(x).  Thus  the  outside  optimizer 
chooses  a  force  level  x  and  an  optimal  and  feasible  use  of  these 
forces  for  each  of  the  objectives,  namely  u1  e  U1(x)  (i  =  1,  ...,  m). 
The  inside  optimizer  makes  his  feasible  choices  v1  e  V1  (i  =  1,  ...,  m) 
to  minimize  the  effects  of  the  outside  optimizer's  choices. 

It  is  shown  in  [3]  that  Problem  (B)  is  equivalent  to  choosing  x 

,1  m  J 

and  u  ,  . . . ,  u  to 


minimize  f(x) 


subject  to 

x  e  X  , 

u1  e  U1(x)  ,  i  =  1,  ...,  m  , 

i.  i  i .  . 

min  g  (x,  u  ,  v  )  >  0  ,  l  =  1,  ...,  m 

•i  „i 
v  eV 

Variants  of  Problem  (B)  analogous  to  those  discussed  above  for 
Problem  (A)  can  also  be  treated. 
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III.  STRATEGIC  OFFENSE  OR  DEFENSE  FORCE  STRUCTURE  OPTIMIZATION 


This  section  presents  an  offense  and  a  defense  force  structure 
model.  Each  is  a  mathematical  program  with  nonlinear  programs  in  the 
constraints.  Both  are  convex  programs.  A  passive  defense  model  is 
presented  in  detail. 

In  the  strategic  offense  model,  the  outside  optimizer  chooses 
minimum-cost  offensive  forces  capable  of  achieving  specified  destruction 
of  various  resources  by  type,  despite  allocation  of  specified  defensive 
forces  of  the  inside  optimizer  to  minimize  the  destruction. 

In  the  strategic  defense  model,  the  outside  optimizer  chooses 
minimum-cost  defensive  forces  capable  of  assuring  specified  survival 
of  various  resources  by  type,  despite  the  optimal  allocation  of 
specified  offensive  forces  of  the  inside  optimizer  to  minimize  the 
surviving  resources. 

Definitions 

Let  the  indexes  of  locations  be  i  =  1,  ...,  p,  of  offensive 
weapon  types  be  j  =  1,  ...,  q,  and  of  defensive  resource  types  be 
k  =  1,  ...,  r.  Let  the  index  of  target  types  be  £  -  1,  ...,  s. 

Define 

y  =  offensive  weapons  of  type  j  targeted  to  location  i, 

z.,  =  defensive  resources  of  type  k  assigned  to  location  i, 

=  parameter  associated  with  destruction  of  target  type  £ 
in  location  i  by  offensive  weapons  of  type  j , 
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0 

^ik  =  Parameter  associated  with  defense  of  target  type  j l  in 
location  i  by  defensive  resources  of  type  k, 

WT  =  value  of  target  type  Si  at  location  i 
f(y)  =  cost  of  providing  offensive  weapons  y  =  ( y _.  ) ( i  =  1, 
•••>  P  >  3  ~  • • • >  Q)  > 


g(z) 


cost  of  providing  defensive  resources  z  =  (z^Xi 
•  •  •  j  p>  k  lj  •  •  •  j  r) . 


1, 


Effectiveness  Function  for  Offense  Optimization 

The  measure  of  effectiveness  for  offense  optimization  is  destroyed 
value  of  type  l  .  One  function  which  is  concave  in  y. .  and  convex  in 


I  '/  exp  ( 
i 


i  ^ik^1 


exp  (  -  E  ai:jyi;j)] 


In  the  strategic  offense  models  the  outside  optimizer  chooses  y^  and 
the  inside  optimizer  chooses  z^,  so  the  function  is  appropriately 
behaved  to  yield  a  convex  program. 


Effectiveness  Function  for  Defense  Optimization 

The  measure  of  effectiveness  for  defense  optimization  is  surviving 
value  of  type  l  .  One  function  which  is  concave  in  z and  convex  in 


y.  .  is 


E  { 1  -  exp  (  -  E  -  exp  (  -  E  a^y^)]} 


In  the  strategic  defense  models  the  outside  optimizer  chooses  z.,  and 
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the  inside  optimizer  chooses  y^,  so  the  function  is  appropriately 
behaved  to  yield  a  convex  program. 

Force  Structure  Optimization  Models 

The  strategic  offense  optimization  model  is  to  choose  y„ 

(i  —  lj  •  •  * > P )  3*  —  •••?  Q)  and  ^  —  •  •  •  >  P?  3  —  1?  •  •  •  >  to 


minimize  f(y) 


subject  to 


min  j  T.  w4  exp(-  £  «-kzik)[l  - 


exp(-  e  a^.y^)] 


s  zik  ±  zk 


>  DX 


where 

D-^  =  specified  destroyed  value  to  be  assured  by  outside  optimizer, 

Z,  =  defensive  forces  of  type  k  available  to  inside  optimizer. 

K. 

The  strategic  defense  optimization  model  is  to  choose  z^ 

(i  =  1,  p;  k  =  1,  ...,  r)  and  y^.C  i  =  1,  p;  j  =  1,  q)  to 

minimize  g(z) 


subject  to 


min  j  S  V/  jl  -  exp(  -  E 


E  y..  <  Y. 
J  ID  -  3 


exp(-  E  a^jy±j )] 


>  S£ 


where 


— 1  =  specified  surviving  value  to  be  assured  by  outside  optimizer, 
=  offensive  forces  of  type  j  available  to  inside  optimizer. 


Y. 

3 
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Properties  of  Models 

The  two  models  are  convex  programming  problems.  They  contain 
convex  programming  problems  in  the  constraints.  Thus  the  procedure 
of  Reference  [2]  can  be  used  for  their  solution. 

Example:  Passive  Defense  Model 

The  outside  optimizer’s  problem  is  to  provide  hardening  and/or 
evacuation  capabilities  for  population  and  hardening  and/or  dispersion 
capabilities  for  industry.  Specified  survival  levels  must  be  achieved 
in  the  face  of  a  nuclear  attack  by  the  inside  optimizer.  These  passive 
defense  measures  are  to  be  supplied  at  minimum  cost  to  protect  against 
attacks  on  either  population  or  industry.  Let  A  -  1  denote  population 
and  A  =  2  denote  industry. 

Define 

=  offensive  weapons  of  type  j  targeted  to  location  i  when 
attack  is  against  target  type  A, 

=  offensive  weapons  of  type  j  available, 

=  hardness  of  target  type  a  at  location  i, 

=  evacuation  capability  of  population  or  dispersal  capability 
for  industry  from  location  i, 

=  scaling  factor  for  damage  by  offensive  weapon  of  type  j 
on  target  type  A  in  location  i, 

=  ‘hardness  parameter  for  target  type  A  at  location  i, 

=  target  evacuation  or  dispersion  parameter  at  location  i, 

=  value  of  the  target  type  A  at  location  i. 


Y. 

3 

4 

4 


a 


Bil 

°i2 

W* 
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J 


Let  the  surviving  value  of  target  type  SL  at  location  i.  for 


offensive  allocations  y^(j  =  1,  q)  be  given  by 


-  exp  (-  8^*^  -  1  -  exp(-  ^  a£,y£,) 


.  SL  SL  SL  9 

The  expression  exp(-  B^z^-  B^2z^2)  9^-ves  t^ie  fraction  of  the  target 
that  is  susceptible  to  attack.  Thus  the  effect  of  increasing  z4^ 
and/or  z^2  is  to  remove  a  portion  of  the  target.  If  z}0  =  0,  then 


i2 


none  of  the  population  is  evacuated  and  if  zf2  =  0,  then  none  of  the 

0 

industry  is  dispersed.  The  variables  z_^  may  be  bounded  from  below, 

_ n 

say  by  to  represent  the  natural  hardness  of  the  population  and 


disperson  of  industry. 
The  expression 


f1  - exp  (-  X  4Ai) 


gives  the  fraction  of 


the  target  which  is  destroyed  by  the  attack  y . .  (j  =  1,  q). 

^  SL  SL 

The  sum  E  a,— y,—  provides  a  measure  of  the  joint  effects 


j=l 


iJ  13 


of  various  types  of  weapons .  The  parameters  scale  weapons  of 
different  yields  to  an  equivalent  number  of  a  standard  weapon. 

A  cost  function  for  hardening  the  targets  might  be 


,  1 
°<zii 


where  c  >  0.  z.,  >  z1, ,  and  d  >  1.  Assume  further  that  the  cost  of 

9  i  1  —  i  r  — 


evacuation  is  linear,  namely,  e^z^2>  ^  similar  cost  function  for 

industrial  protection  is  assumed. 
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- ^ 


The  overall  problem  of  providing  defensive  forces  at  minimum  cost  to 

achieve  weighted  population  survival  r^  and  weighted  industry  survival  r^ 

1  1  2 
is  to  choose  z±1  (i  =  1,  p),  zH9  (i  =  1,  p),  z_.  (i  =  1,  p), 


Ji2 


il 


i 9  _  •••>  p)>  yjj  ^  1»  •••»  p;  3  1»  •••»  y±j  c1  -  i>  •••>  p 


i2 

3=1,  q)  to 

minimize  E  cA(2t.  -  zt.  )Vi  +  E  e?z?9 

i=l  ax  ax  i=i  i 

P  o  2  _o  d^  p  2  2 
+  E  c  (Zil  -  zA1)  +  E  eizi2 


p  1/_1  -1  vd1  p 


1  1 


i=l 


i=l 


subject  to 


minimum  E  W^. 

y1.  i=1 

13 


i  - 


/  ol  1 
exp  (-  3ilzil 


-  -  exp  (-  E  cu.yK)! 

pi2  i2'[_  F  v  j=1  13  13  J 


P  i 

E  y. .  <  Y.  , 
i=l  *3  ~  3 


3  _  i,  •  •  • ,  q 


>  r' 
-  1 


minimum  E  W? 
2  i=l  1 

yij 


1  - 


,  02  2 
exp  (-  P^z^ 


o  2  T  q  2  2  1 

"  3i2Zi2)[1  ‘  eXP  (~  aijyij 


E  y  t .  <  Y .  ,  3  =  1,  •  •  • ,  q 


i=l 


13-3 


—  r2 


zil  —  zil  ’  i  1>  • • • >  P 


2-2 
al  —  il  * 


i  —  lj  •  •  • ,  P 
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IV.  STRATEGIC  BOMBER  FORCE  STRUCTURE  AND  BASING  OPTIMIZATION 

This  model  is  an  extension  of  the  model  presented  in  Reference  T4] . 
The  problem  is  to  provide  a  minimum-cost  force  structure  and  basing 
of  strategic  bombers  which  will  enable  the  outside  optimizer  to 
achieve  specified  survivability  in  the  face  of  a  submarine -launched 
ballistic  missile  (SLBM)  threat.  The  outside  optimizer  chooses  bomber 
forces  and  basing,  and  the  inside  optimizer  chooses  allocations  of 
submarines  to  launch  areas  and  targeting  of  SLBMs  to  bases  to  minimize 
surviving  bombers . 

Let  i  =  1,  p  index  the  missile  launch  area,  j  =1,  . ..,  q 

index  the  bomber  bases  and  t  =  0,  1,  ...,  r  index  time.  Define 


x. 

3 


f(xx. 


•>  v 

gt(x.) 

3 


Mi 


M 


number  of  bombers  assigned  to  base  j  , 

cost  of  providing  bombers  x1 ,  . . . ,  x  , 

number  of  bombers  on  base  j  at  time  t  when 

the  number  at  t  =  0  is  x.  , 

3 

number  of  missiles  launched  from  area  i  which 
are  targeted  to  arrive  at  base  j  at  time  t  , 
probability  that  a  missile  launched  from 
area  i  hits  base  j  at  time  t  , 
missile  flight  time  from  launch  area  i  to 
base  j  , 

number  of  missiles  at  launch  area  i  , 
total  number  of  missiles  available. 
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If  the  objective  of  the  attacker  is  to  minimize  the  number  of  bombers 


surviving,  his  problem  is  to  choose  v5.  and  to  minimize 


q 

Z 

j=l  t=0 


t 

v .  . 

r  p  ij 

E  gj( xd)  *  (i  -  p*j> 


C-1  s 

TT  (1  ~ 


13 


s=0 


13 


subject  to 

q  t+t 

E  vu  <  Mj  ,  1  =  1 . P,  t  =  0,  ....  r 

3=1 

P 

Z  M.  <  M 
i=l 


Let  S  be  the  required  number  of  surviving  bombers  after  absorbing 
the  optimal  attack.  The  overall  problem  is  to  choose  x^,  , 

vj.  (i  =  1,  p;  j  =  1,  . q;  t  =  0,  r),  and  (i  =  1,  p.) 

to 

minimize  f(x^,  x  ) 


subject  to 


P  v..  t-J 

_  t,  X  /T  t  X  13 

Z  g,.(X^)  TT  (l  -  Pij) 


q  r 

minimum  Z  _  _  „ 

t  M  3=1  t=0  3  3  i=l 

v  .  . ,  M.  J 
13’  1 


tt  d-P-j)^ 
s=0  J 


q  t+t.. 

E  v .  .  J  <  M  •  ,  i  =  1,  .  . . ,  p,  t  =  0,  . 
3=1  3 


P 

EM.  <  M 
i=l 


J 


>  S 
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The  optimization  problem  in  the  constraint  is  a  convex  nonseparable 
nonlinear  program.  A  detailed  example  of  a  special  case  of  the  inside 
mathematical  program  is  given  in  Reference  [4] .  If  the  functions 

are  concave  in  x^ ,  then  the  constraints  satisfy  the  conditions 
necessary  for  a  convex  program.  Finally,  if  f(x^,  ...,  x  )  is  convex, 
the  overall  problem  is  a  convex  program. 
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V.  STRATEGIC  DEFENSE  OPTIMIZATION  TO  ACHIEVE 

SPECIFIED  POST-ATTACK  PRODUCTION  CAPABILITIES 

This  model  is  described  in  more  detail  in  Reference  [5].  It  is 
summarized  here. 

One  of  the  planning  considerations  for  defense  against  a  nuclear 
attack  is  the  assurance  of  sufficient  surviving  economic  capacity  to 
support  the  surviving  population.  The  present  model  contains  an 
aggregate  representation  of  the  nation-wide  production  base.  The 
effect  of  strategic  defenses,  consisting  of  both  active  defenses  (for 
example,  anti-ballistic  missiles)  and  passive  defenses  (for  example, 
population  shelters),  against  offensive  weapons  is  modeled.  The 
defender  chooses  a  minimum-cost  mix  of  active  and  passive  components 
which  ensures  that  specified  post-attack  production  capacities  will 
survive  after  an  optimized  attack. 

The  country  is  partitioned  into  j  =1,  ...,  n  geographic  regions, 
allowing  consideration  of  varying  population  and  production-base 
densities.  The  general  model  allows  i  =  1,  ...,  m  different  economic 
sectors  in  each  geographic  region.  Each  economic  sector  is  charac¬ 
terized  by  a  Cobb-Douglas  production  function.  Also,  k  =  1,  ...,  p 
denotes  the  different  types  of  defensive  resources  and  1  =  1,  . . . ,  Q 
the  different  types  of  offensive  weapons. 

Define 

x  .  =  number  of  defensive  resources  of  type  k  assigned 
to  region  j  , 
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=  number  of  offensive  resources  of  type  £  targeted 
on  region  j  in  an  attack  on  economic  sector  i  , 

=  number  of  offensive  weapons  of  type  £  . 


The  post-attack  production  function  (in  terms  of  value  added)  in 
economic  sector  i  in  region  j  is  assumed  to  be 


H.  . 
31 


-  Z  a"  x., 
k--*l  ljk  :,k 


•v  b".y. 

fiy  J £ 


1  -  e 


]  -  c 


CL  *  * 


iaii 


-  £  a*  x.. 
k:-l  Jjk  Jk 


?  1 K  i 

■  T.  b  .  .  .  v  . 
1J£ 


1  -  c 


1  -  c 


* 


1  I 

.  J]  a  ..  x., 

k-i Jl; 


£  ,  J.  i 

■  •  h  .  .  .V  .  . 

£«1  ^  ^ 


1  -  e 


1  -  c 


8  .  . 
3.1 


If  v*  =  0,  £  =  1,  . . . ,  q,  then  the  production  function  has  the 
J  £ 

standard  Cobb-Douglas  form 


a-  •  6. . 

H..  K.1J  L'V 

33  33  13 


where  H_._.  represents  the  technological  efficiency,  is  the  capital 

base,  and  L. .  is  the  labor  base.  The  exponents  a- •  And  8. .  are  the 
i 3  33  33 

elasticities  of  value  added  with  respect  to  capita]  and  labor, 

respectively. 
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The  expressions  of  the  form 


1  - 


-  E  a .  x., 

k_,  13k  3* 

e  k-J-  \  1 


-  e 


modify  the  efficiency,  the  capital  base,  and  the  labor  base  as  a 
function  of  the  offense  and  defense  allocations.  The  expression  is 
convex  in  v^,  concave  in  x^,  and  has  the  asymptotic  properties 
expected  for  the  physical  processes  being  modeled.  The  parameters 

a..,  and  b..„  can  be  estimated  from  detailed  analyses, 
lgk  lji 

The  aggregate  production  for  the  whole  country  is  the  sum  of  the 
production  across  geographic  regions. 


P  n 

The  cost  of  defensive  resources  is  taken  to  be  E  E  c.,x.,  > 

k=l  j=l  Jk  Jk 

where  c.,  is  the  unit  cost  of  defensive  resource  k  in  region  j  . 

3k 

The  requirements  for  surviving  post-attack  production  capacity 
are  given  by  r.(  i  =  1,  ...,  m). 

The  overall  model  is  to  choose  x.^  (  j  =  1,  ...,  n;  k  =  1,  ...,  p) 

a  nd  v .  .  (  l-l,  ...}  m ;  j  1,  ...,  n,  1,  ...,  q  )  to 

3* 

P  n 

minimize  E  £  c.,x 
k=i  j=i  ]k 


*> 


17 


subject  to 


n 


minimum  EH.. 

i  j=l  *3 
v.  „  J 


-  Z  a »  x*vf  -  i  bH.  ,v* 

!_.ek=1  3kl  Aal  13£31 


•  k!^ 
13 


8,, 


n 


£  -X 

.  k=i  i3k  3k  . 

1  -  e  \1  -  q 


•  q 

£=rijrje 


-  E  b'f.  v1 


ID 


-  E  a^  x  /  -  y  bi*.  v4 

X  -  .  k=k  13k  3k(l  -  e  ^  131  3£ 


E  Vi  (  -  V(.  ’ 

j=l  3t  1 


^  i>  •••>  q 


>  ri  ,  i=l, . . .  ,m. 


The  overall  model  is  a  convex  program.  The  inside  optimization  problems 
are  convex  nonseparable  nonlinear  programs . 
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VI.  OPTIMIZATION  OF  WEAPON  MIX  AND 

TARGETING  FOR  ATTRITION  PROCESSES 


A  conventional  warfare  model  is  formulated  for  the  attrition  of 
a  heterogeneous  mix  of  weapons.  The  problem  for  the  outside  optimizer 
is  to  choose  a  minimum-cost  mix  of  weapons  and  their  optimal  targeting 
patterns  against  the  optimal  targeting  of  the  weapons  of  the  inside 
optimizer.  The  force  size  must  be  sufficient  to  allow  the  outside 
optimizer  to  achieve  a  specified  differential  value  of  weapons  at 
the  end  of  the  engagement. 

Let  the  index  i  =  1,  . . . ,  m  denote  weapon  type  for  the  outside 
optimizer  and  j  =1,  . . . ,  n  denote  weapon  type  for  the  inside 
optimizer.  Define 

x^  =  initial  number  of  outside  optimizer’s  weapons  of  type  i  , 

y_.  =  initial  number  of  inside  optimizer's  weapons  of  type  j  , 

u^j  =  number  of  outside  optimizer's  weapons  of  type  i  assigned 

to  fire  on  inside  optimizer's  weapons  of  type  j  , 

.  v„  =  number  of  inside  optimizer's  weapons  of  type  j  assigned 
to  fire  on  outside  optimizer's  weapons  of  type  i  , 
b^.  =  effectiveness  of  outside  optimizer's  weapon  type  i  on 
inside  optimizer's  weapon  type  j  , 
r_  =  effectiveness  of  inside  optimizer's  weapon  type  j  on  outside 
optimizer's  weapon  type  i  . 

In  classical  Lanchester  theory  the  instantaneous  attrition  to 
a  weapon  by  type  i  from  an  opposing  weapon  of  type  j  is  given  by 
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x.  =  -  r . .y . 
i  ]i  ] 


or 


x .  =  -  r  . .  x .  y .  , 
i  ji  a/j  ’ 

the  first  equation  being  the  "square"  law  and  the  second  the  "linear" 
law.  More  generally, 

•  2-p 

x.  =  -  r . .x .  *  y.  , 
i  31  i  j 

where  1  <  p  <  2,  gives  the  "p"-law,  with  the  square  and  linear  laws 
being  special  cases. 

t 

The  heterogeneous  Lanchester  attrition  considers  more  than  one 
weapon  type  for  each  opponent.  The  heterogeneous  analog  to  the 
p-law  is 

2-Pi  n 


xi  =  -  x 


E  r  .  .  v  . .  , 
1  •  1  ji  ]i  ’ 

1=1  J  J 


where  1  <  p.  <  2,  i  =  1,  . . . ,  m. 


In  conventional  Lanchester  analysis  the  parameters  b^  (and  rj^) 
are  instantaneous  attrition  rates.  For  purposes  of  the  model  here  it 
is  useful  to  consider  these  parameters  to  be  the  kill  rates  over  some 
positive  length  of  time.  Thus,  b^.  might  be  the  average  effectiveness 
of  a  weapon  type  i  firing  at  a  weapon  type  j  for  a  day,  a  week,  or 
longer.  The  number  of  surviving  weapons  of  type  i  for  the  outside 
optimizer  at  the  end  of  the  length  of  time  is 


x. 

1 


2-Pi 


n 

E 

j=l 


r .  .v. . 
li  11 
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Similarly  the  surviving  weapons  of  type  j  for  the  inside  optimizer 


are 


m 


£ 

i=l 


b. .u. . 
il  1 1 


3 


where  1  <  <  2,  j  =  1,  . m. 

If  0^  and  p_.  are  the  values  of  the  surviving  weapons  for  the 
outside  optimizer  and  for  the  inside  optimizer,  respectively,  then 
the  model  for  finding  minimum-cost  forces  subject  to  the  difference 
in  value  of  surviving  weapons  being  at  least  r  is  to  choose 
x  =  (Xp  xm),  u  =  (u^. ),  and  v  =  Cvji)(l  =  1,  •  ••>  m; 

j  =  1,  . . . ,  n)  to 


minimize  f(x) 
x 


subject  to 


1“ 

[ m  r 

2-p  •  n 

n  r 

2-q-  m  f 

max 

[ueU( x) 

min]  £  0 . lx. 

veV (i=l 

-  x.  £  r.  .v. . 

i  j=1  H  1  ij 

-  £  Pj  |y 

j=i  1 L 

.  -  y .  E  b..u.  .1} 

1  yl  i=1  H  iljj 

u 

A| 

where 


n 

U(  x)  —  ]  ^j  >  0:  .  £  Uij  x^ ,  i  1  j  •••?  m 


m 


V  =  ]  v . .  >0:  £  v..  <  y. ,  j  =1,  n>  , 


li  - 


i=l 


li  ~  1 


and  where  f(x)  is  the  cost  of  forces  x  . 


21 


Applying  the  result  in  Reference  [3],  the  above  optimization 
model  is  equivalent  to  the  following: 

minimize  f(x) 
x,  u 

subject  to 

n 

E  u. .  <  x.  ,  i  =  1,  . . . ,  m 

.  ,  lj  —  i  ’  ’  5 

3=1  J 


u^  >0  ,  i  =  1,  m;  j  =  1,  n 


The  inside  mathematical  program  is  a  linear  program.  Further,  the 
inside  objective  function  is  concave  in  (x,u).  Thus,  if  f(x)  is 
convex,  the  program  is  convex.  There  is  a  possibility  that  an 
'’overkill”  of  forces  will  occur  if  the  attrition  parameters  are  not 
small  relative  to  the  number  of  weapons.  The  numerical  solution 
of  the  problem  should  be  examined  for  this  eventuality. 
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VII.  OPTIMIZATION  OF  AIRCRAFT  DEPLOYMENT  AND  SORTIE  ALLOCATION 

The  purpose  of  this  model  is  to  choose  optimal  Blue  aircraft 
deployments  to  theater,  and  Blue  and  Red  optimal  sortie  allocations, 
assuring  that  Blue  achieves  specified  differences  of  cumulative  ground 
and  air  firepower  at  the  end  of  a  war.  The  war  consists  of  two  periods. 
Red  aircraft  deployments  are  fixed.  The  sortie  allocations  are  to 
combat  air  support  (CAS)  and  air-base  attack  (ABA). 

Let  t  =  1,  2  index  time,  and 

x^  =  number  of  Blue  aircraft  deployed  to  theater  for  use  in 
period  t 

yt  =  number  of  Red  aircraft  deployed  to  theater  for  use  in 
period  t 

ul,u2  =  num^er  Blue  aircraft  assigned  to  CAS  and  ABA, 
respectively,  in  period  1 

v^,v^  =  number  of  Red  aircraft  assigned  to  CAS  and  ABA, 
respectively,  in  period  1 

a  =  kill  parameter  for  Blue  ABA  killing  Red  aircraft 
B  =  kill  parameter  for  Red  ABA  killing  Blue  aircraft 
f  =  firepower  per  Blue  CAS  sortie 
g  =  firepower  per  Red  CAS  sortie 
c.  =  cost  per  Blue  aircraft  introduced  at  time  t 
S  =  required  cumulative  firepower  difference  by  time  t. 

Destruction  of  aircraft  is  a  function  of  the  number  of  opposing 
ABA  aircraft.  The  attrition  of  Blue  aircraft  is  given  by 

-Bv2 

x(l  -  e  )  , 
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and  similarly  for  Red.  Both  sides  allocate  all  of  their  aircraft 
to  CAS  in  the  second  period,  since  aircraft  allocated  to  ABA  would 
have  no  effect  due  to  the  war  having  ended. 

The  two-period  optimization  model  is  to  choose  x^,  x2,  u^,  u2, 
v2  to 

minimize  c^x^  +  C2X2 

Xl>  x2 

subject  to 

max 

ul+u2<xi 

Applying  the  result  of  Reference  [3]  the  problem  can  be  written  as 
to  choose  Xp  x2,  u  ,  u2,  v.^  v2  to 


min  +  (x]_+x2-x]_(l  ~  e 


-9v, 


)] 


-au?  , 

-  g[v1  +  (y1+y2-y1(l  -  e  ))]| 


>  S2 


minimize  c..x^  +  c2x2 
x1,x2,u1,u2 


subject  to 


mm 

vl’v2 


-Svr 


(  2 
^f[ux  +  (x1+x2-x1(l  -  e  )] 


-au2  | 

sO-l  +  (y1+y2-y1(1  -  e  )]> 


vi  +  v2  ^  yi 
U1  +  u2  —  xl 


>  S2 


The  overall  program  is  convex,  and  the  inside  program  is  also  convex. 
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